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We consider a financial contract that delivers a single cash flow given by the termi- 
nal value of a cumulative gains process. The problem of modelling and pricing such 
an asset and associated derivatives is important, for example, in the determination of 
optimal insurance claims reserve policies, and in the pricing of reinsurance contracts. 
In the insurance setting, the aggregate claims play the role of the cumulative gains, 
and the terminal cash flow represents the totality of the claims payable for the given 
accounting period. A similar example arises when we consider the accumulation of 
losses in a credit portfolio, and value a contract that pays an amount equal to the 
totality of the losses over a given time interval. An expression for the value process 
of such an asset is derived as follows. We fix a probability space together with a 
pricing measure, and model the terminal cash flow by a random variable; next, we 
model the cumulative gains process by the product of the terminal cash flow and an 
independent gamma bridge process; finally, we take the filtration to be that gener- 
ated by the cumulative gains process. An explicit expression for the value process is 
obtained by taking the discounted expectation of the future cash flow, conditional 
on the relevant market information. The price of an Arrow-Debreu security on the 
cumulative gains process is determined, and is used to obtain a closed-form expres- 
sion for the price of a European-style option on the value of the asset at the given 
intermediate time. The results obtained make use of various remarkable properties 
of the gamma bridge process, and are applicable to a wide variety of financial prod- 
ucts based on cumulative gains processes such as aggregate claims, credit portfolio 
losses, defined-benefit pension schemes, emissions, and rainfall. 
Key words: Asset pricing; insurance claims reserves; credit portfolio risk; cumula- 
tive gains, gamma bridge process; beta distribution; option pricing; reinsurance. 
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I. INTRODUCTION 



There are a number of problems in finance and insurance that involve the analysis of 
accumulation processes — that is to say, processes representing cumulative gains or losses. 
The typical setup is as follows. We fix an accounting period [0, T] , where time denotes the 
present. At time T a contract pays a random cash flow X T , which is assumed to be positive 
and given by the terminal value of a process of accumulation. In the case of an insurance 
contract, for example, we consider the situation where a number of claims are made over the 
accounting period, and are then paid at T. The random variable Xt represents the totality 
of the payments made at T in settlement of claims arising over the accounting period. The 
problem facing the insurance firm is the valuation of the random cash flow. Let us write 
{St} for the value process of the contract that pays Xt at T, and {Tt\ for the filtration 
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representing the flow of information available to market participants, and Q for the pricing 
measure, which we assume to have been established by the market. Then the value at t of 
the contract that pays X T at T is 

S t = P tT E[X T \f t ], (1) 

where E[— ] denotes expectation with respect to Q, and P tT denotes the discount factor, 
which for simplicity we take to be deterministic. One can interpret S t as the reserve that 
the insurance firm requires at t to ensure that Xt will be payable at T. Alternatively, 
one can view S t as the amount that would have to be paid at t in order for the insurance 
firm to relieve itself of the obligation to pay Xt, that is to say, to commute the relevant 
claims. Similarly, the cost C tT at t of a simple stop-loss reinsurance contract that pays out 
(X T — K) + at T for some fixed threshold K is given by 

C tT = P tT E[(X T - K) + \T t \. (2) 

We shall assume that {Ft} is generated by an aggregate claims process {£t}, where for 
each t the random variable & represents the totality of claims known at t to be payable at T. 
The problem can then be stated as follows: given the history of claims over the accounting 
period up to time t, what is the appropriate reserve to allocate for settlement of these and 
any future claims arising in the accounting period? To obtain a solution to the problem we 
need to specify the aggregate claims process, then work out the reserve process {S t }. Once 
we have the reserve process, we can value various types of reinsurance contracts. 

Another example of an accumulation process comes from credit risk management. We 
consider a large credit portfolio, and let Xt denote the value of the accumulated losses 
at T. For instance, at time a credit-card firm has a large number of customers, each 
with an outstanding balance payable in the accounting period. If a customer does not pay 
the balance by the required date, they will be deemed to be in default, and a loss will be 
registered. The random variable X T will denote the totality of such losses. We assume 
that once a customer is in default, no further payments are made by that customer (this 
assumption can be relaxed in a more sophisticated model). The problem facing the credit- 
card firm is to determine what reserve policy to maintain, and what premium to charge over 
the base interest rate, to ensure that funds will be in hand to cover the default losses. 

The purpose of this paper is to present a modelling framework for accumulation processes, 
and to establish explicit formulae for the associated valuation processes. In particular, we 
shall assume that {£ t } takes the form 

& = X T itT, (3) 

where {•jtr} is a gamma bridge over the interval [0,T], independent of Xt- The motivation 
for the form of the accumulation process indicated above arises in two distinct lines of 
enquiry. The first relates to the idea that the gamma process might be used as a basis 
for describing the aggregate losses associated with insurance claims. This idea dates to the 
work of Hammersley (1955), Moran (1956), Gani (1957), Kendall (1957), and others, in 
connection with the theory of storage and dams. Moran (1956), in particular, observed that 
the amount of rainfall accumulating in a dam can be modelled by a gamma process, and 
Gani (1957) pointed out the relevance to insurance, the argument being that providing that 
the portfolio of events insured is sufficiently large, one can think of the arrival of claims 
as being analogous to the accumulation of dam rain. The gamma process has since then 
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been investigated by Dufresne et al. (1991), Dufresne (1998), Dicksen & Waters (1993), and 
others, as a model for aggregate claims. 

Let us therefore consider what results if we model the aggregate claims process as a Q- 
gamma process. In other words, suppose we set £ t = Kj t , where k is a constant and {jt} is a 
standard gamma process under Q, with mean and variance mt (see Section [IT] for definitions). 
It follows that = Xr7tT where Xt = kjt and the process {~itr} defined by •ytr = It/ It 
is a standard gamma bridge over [0, T\. Moreover, by virtue of the special properties of the 
gamma process, we find that Xt is independent of {7<t}- We see that in the Q-gamma 
model the aggregate claims process is the product of a gamma-distributed terminal cash 
flow and an independent gamma bridge. One can think of the gamma bridge as representing 
that aspect of the aggregate claims process that has no bearing on the terminal result. We 
are thus led to a multiplicative decomposition of the accumulation process into the product 
of a "signal" Xt and an independent "noise" {7<t} carrying no information about Xt- 

For such processes we are able to apply the techniques of information-based asset pricing 
developed in Brody et al. 2007a, 2007b, Hughston & Macrina 2007, Macrina 2006; and 
Rutkowski & Yu 2007. Indeed, through this second line of enquiry one is led to consider the 
more general situation where the terminal cash flow, instead of being gamma distributed, 
has a generic a priori distribution, and the claims process takes the form ([3]). The additive 
decomposition of the market information process in the case of the Brownian bridge noise 
considered in the references cited above is natural from the viewpoint of nonlinear filtering 
theory. The product representation of the gamma information process is equally natural, 
since many properties of the Brownian bridge that hold additively have striking multiplica- 
tive analogues for gamma bridges (Emery & Yor 2004, Yor 2007). The resulting model 
for the aggregate claims process is remarkably tractable, and we are able to derive explicit 
formulae both for the claims reserve process, and for the valuation of reinsurance contracts. 

The paper is organised as follows. In Sections II, III, and IV, we outline a number of the 
properties of gamma processes and gamma bridges. The material covered in these sections 
is for the most part well known. However, since it is not easy to locate a systematic but 
elementary treatment of the gamma process and the associated bridge process, it will be 
useful to present some of the details here for the benefit of general readers. At the same 
time, we establish our notation and some results that will be applied in later sections. In 
Section |V] we derive an explicit expression for the value process of a contract that delivers 
the cash flow Xt at time T, when the market filtration is generated by the accumulation 
process ([3]). We show in Proposition @] that {^} has the Markov property, and then use the 
Bayes theorem to determine the conditional density of Xt-, and finally the value process, 
which is given in Proposition [51 By use of the conditional density we are also able to obtain 
an expression for the value process of a simple stop-loss reinsurance contract. In Section IVT1 
we consider the valuation of general reinsurance contracts. In particular, we derive a formula 
for the value at time of a contract that at some fixed time t gives the contract holder the 
option to commute the claim Xt by paying a fixed amount K at t. Such a contract takes the 
form of a European call option on the value of the reserve at t. An Arrow-Debrue method 
is introduced to simplify the calculations. The resulting formula for the option value is 
expressed in terms of the cumulative beta distribution. We examine in Section IVHI the case 
where Xt takes discrete values. When Xt is a binary random variable, the problem of option 
pricing can be solved completely. In Section IVIIII the material of Section IVII is extended to 
determine an expression for the price process of an option on the value of an aggregate claim. 
In Section IIXI we conclude by returning to the case where Xt has a Q-gamma distribution. 
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II. GAMMA PROCESSES AND ASSOCIATED MARTINGALES 

We fix a probability space (Q, J 7 , Q). In our applications Q will always denote the pricing 
(risk-neutral) measure, but the material in this section, and the following two, does not 
depend on this interpretation. Equalities and inequalities among random variables are to 
be understood as holding except possibly on sets of measure zero. By a standard gamma 
process {-ft}o<t<oo on (f2,jF, Q) with growth rate m we mean a process with independent 
increments such that 70 = and such that the random variable j t has a gamma distribution 
with mean and variance mt. More precisely, writing G(x) = Q[ 7t < x ] f° r the distribution 
of 7t, and writing g(x) = dG(x)/dx, we have 

g(x) = i {x>0} —^- (4) 

for the density of j t . Here T[a] is the standard gamma function, which for a > has the 
Eulerian representation: 

POO 

T[a] = / x a - l e~ x dx. (5) 
Jo 

It follows from the identity F[a + 1] = aT[a] satisfied by the gamma function that 

E[ 7t ] = mt, (6) 

which justifies the interpretation of the parameter m as the mean growth rate of the process. 

A straightforward calculation shows that the characteristic function for the gamma pro- 
cess is given by 

valid for t > and for A G C such that Im(A) > —1, from which the higher moments of ^ t 
can be deduced. We note that E [7 t 2 ] = mt + m 2 t 2 , and hence that Var[7 t ] = mt. It follows 
as a consequence of the independent increments property that Cov[7 t , 7J = mt for u > t. 

An alternative expression for the characteristic function is given by the Levy-Khinchine 
representation E [e lA7t ] = e~*^ A ) for Im(A) > — 1, where 

POO 

jjj(X) = m ln(l - iA) = / mx~ X Q~ x (l - e iXx ) dx, (8) 

Jo 

which shows that the Levy density associated with the gamma process is given by mx~ x e~ x 
for x > (see, e.g., Protter 2005). 

By use of the independent increments property we deduce that for u > t > and for 
a, b E C with Im(a + b) > — 1 and Im(fe) > —1 we have: 

= E [ e i(a+b)7t ] E [e ib(7u " 7t) ] 

= I I (9) 

[1 -\{a + b)] mt (1 -i6)m(u-t) - 1 > 
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In particular if we set —a = b = A, we see that ^ u — 7$ is gamma-distributed with parameter 
m(u — t). It follows that the increments of {7 t } have a time-homogeneous probability law 
in the sense that j u +h — Jt+h has the same distribution as 7„ — j t - 

Using the independent increments property it is a straightforward exercise to deduce that 
the processes {*yt — Tnt} and {^ 2 — 2mt r )t + mtimt — 1)} are martingales. More generally, 
for a > — 1 the process {L t } defined by 

L t = (1 + a) mt e- ait (10) 

is a martingale, which can be verified by use of Q- We refer to this process as the exponential 
gamma martingale. It follows, by consideration of the corresponding power series in a, 
that for each term in the series we are able to obtain a martingale involving a polynomial 
expression in the gamma process. Suppose for n G IN and k G R we define the so-called 
associated Laguerre polynomials {Ln\z)} by setting 

rT 

LW(z)=z- k e*^(z n+k e-*). (11) 

Thus, we have L {k \z) = -z + k + 1, L (k \z) = \[z 2 - 2{k + 2)z + (k + l)(fc + 2)], and so on. 

The standard Laguerre polynomials, given by L n (z) = Ln\z), have the property that if Z is 
a standard exponentially distributed random variable, then E [L n (Z)L n >(Z)] = for n ^ n' 
(cf. Wiener 1949). More generally, if Z has a gamma distribution with parameter k + 1, 
i.e. such that Q[Z < z] = /* x k e~ x dx/T[k + 1], for k > -1, then E[L {k) (Z)L {k) (Z)] = 
for n 7^ n' . The significance of the associated Laguerre polynomials in the present context 
arises from the identity 

oo 

(l + a) h e- a = j^Lf- B )(z)a» ) (12) 

n=0 

valid for |a| < 1 and h > (Erdelyi 1953), which gives us the required series expansion of the 
exponential gamma martingale in powers of a. In particular, by setting h = mt and z = % 
in equation (1T2]) . we are able to deduce that for each value of n the process {Li m * (7*)} is 
a martingale (cf. Schoutens 2000). For example, we have 

Lt^ilt) = -(it - mt), 4 m '" 2) (7t) = \[l 2 t ~ 2mt lt + mt(mt - 1)]. (13) 

So far we have confined the discussion to the case of the "standard" gamma process, 
for which Epy$] = mt and Var[7 t ] = mt, for some value of m. We note that the ratio 
(E[7i]) 2 /Var[7 t ] is dimensionless, and hence that m has the units of inverse time. For any 
fixed m we can choose the units of time so that m — 1 in those units (this is done implicitly, 
for example, in Yor 2007). We shall, however, take the units of time as fixed, and m as a 
model parameter. 

For many applications it is useful also to consider a broader family of gamma processes, 
labelled by two parameters, which we shall call "scaled" gamma processes. By a scaled 
gamma process with growth rate /1 and spread a we mean a process {T t } < t<oo with inde- 
pendent increments such that r = and such that T t has a gamma distribution with mean 
fit and variance a 2 t, where fi and a are parameters. Defining m = f 2 /^ 2 and k = cr 2 /f-, we 
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have fi = Km and a 2 = K 2 m. One can think of m as a "standardised" growth rate, and k as 
a "scale" . The density of T t is then given by 

g Vt {x) = 1 {X>0} ^—^ . (14) 

It is straightforward to check that if {T t } is a scaled gamma process with standardised 
growth rate m and scale k, then {/t _1 r t } is a standard gamma process, with growth rate m. 

Now suppose that {7 t } is a standard gamma process on (ft, J 7 , Q), let {Q t } denote the 
filtration generated by {jt}, and let Q* denote the measure on (fl, Qt), for some fixed T, 
defined by the likelihood ratio 



*r mr exp(— 7T ) (15) 



T 



K 



for some k > 0. Then {7t}o<t<T is a scaled gamma process on (Q, Qt,Q*), with scale 
parameter k. Thus, E[7 t ] = mi, Var[7 4 ] = mi, E*[ 7 t] = unit, and Var*[7 t ] = K 2 mt. This can 
be established by working out the joint characteristic function under Q* of the increments 
It ~ 7s, 7s ~ 7sn 7si ~ 7s 2 , • • • , 7s n -i ~ 7s n for T > t > s > s 1 > s 2 > ■ ■ ■ > s n for 
each n G IN, and showing that it factorises. We note that the change-of-measure density 
martingale arising in this example is obtained by taking the standard gamma exponential 
martingale (fit)]) defined above, and setting a = (1 — k)/k. 

The gamma process has been used as the basis of a number of different asset pricing 
models; see, for example, Madan & Seneta (1990), Madan and Milne (1991), Heston (1995), 
Madan et al. (1998), Carr et al. (2002), and Baxter (2007). 



III. GAMMA BRIDGE PROCESSES 

Let {7t}o<t<oo be a standard gamma process with growth rate m, and for fixed T define 
the process {7 tr }o<t<T by setting 

7tT = ^- (16) 

7t 

Then clearly 7 t = and ^tt — 1- We refer to {7tr}, thus defined, as the standard gamma 
bridge over [0, T] associated with the gamma process {7t}- More generally, we refer to any 
process having the law of {7^} as a standard gamma bridge over [0, T\. It can be shown 
that the random variable ^tT has a beta distribution. In particular, we have the following: 

Proposition 1 The density function of the random variable is given by 

y-mt-l^ _ y\m(T-t)-l 

M = WD B[mt, m(T — t)] ' 



where 



_ T[a}T[b} 

5M " r>T&T (18) 
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Proof. First we note that 



It . 
— < V 
It 



it < y 



It ~lt l-y 



(19) 



Since j t and 7t — 7t are independent, and 74 has a gamma distribution with parameter mt, 
we have 



It < y 



It ~lt 1 - y 



E 



It < T^—(lT-lt) 

l ~y 



It < YZ~^ (Tt - 7t) 



7t - 7t 



T[mt] 



E 



i-r^(7T-7t) 



(20) 



Therefore, the corresponding density is given by 



7t . 

— < y 

7t 



1 



1 



1 



1 



{0<»<1} 



{0<y<l} 



T[mt] 
1 

T[mt] 



E 



E 



d rj^iiT-it) 



x mt-l e -x dx 



dy 

It ~lt I y 



[l-yY \l-y 



rnt—l 



{0<j/<1}' 



T[mt] 



■mt— 1 



E 



(21) 



Now, since 77- — 74 has a gamma distribution with parameter m(T — t), for the expectation 
appearing in the line just above we obtain 



E 



(7T-70 m 'e"^ {7T " 7t) 



T[m(T-t)} J 

1 r°° 

T[m(T - 1)) I 
(1 - y) mT rco 



T[m(T - 1)} J 
T[mT] 



x^-^-^dx 
u mT - l e- u du 



r[m(T - 1)} 



(i-y) 



mT 



(22) 



where in the last two steps we make the substitution x = u(l — y) and use formula (jHJ). 
Putting this result back into (12~II) . we obtain ( TT71) . as desired. □ 



Let us calculate the moments of 7^. Bearing in mind the integral representation 



B[a,6] 



./o 



(23) 
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for the beta function, we deduce that 

E[ ^ t]= B[mtMT-t)] (24) 

for n > 0. By use of ({TBI along with the identity F[a + 1] = aT[a] we find that IE^y] = t/T 
and that IE [7^] = t(mt + l)/T(mT + 1). It follows in particular that 

Var[ 7tT ] = J} T - %- . (25) 
L/ J T 2 (l + mT) y ' 

It is interesting to observe that the expectation of -y t T does not depend on the growth rate 
m, and that the variance of jtr decreases in increasing m. 

More generally, let us define the Pochhammer symbol by writing (a)o = 1 and (a)k = 
a(a + l)(a + 2) • • • (a + k — 1). Then we find that the moments of jtT are given by the 
expression Epy^p] = (mt) n /(mT) n , and for the corresponding central moments we obtain 

E[( ltT -E[ ltT ]) n ] = (-1) F(-n,mt;mT;T/t), (26) 



where F (a, b; c; z) = Yl'k=o{ a )kip)kZ k /[k\(c)k] is the hypergeometric function (Erdelyi 1953). 



IV. FURTHER PROPERTIES OF GAMMA BRIDGES 



It is a remarkable property of the gamma process and the associated gamma bridge that 
the processes { r y u } u >T and {7tT}o<t<T ar e independent. In particular, the random variables 
7t and ^ t T = It I It are independent for < t < T. This property allows us to verify 
straightforwardly that {74} has the Markov property. To show that {74} has the Markov 
property we need to verify for a > that 

(27) 



Q [it < a>\ls, 7si, 7s 2 , • • • , 7sJ = Qbt < a\j s ] 
for all t > s > Si > S2 > ■ ■ ■ > s n > 0, and for all n > 1. But clearly, 



[it < a|7->7«i>7* a >--->7*n] 



It < a 
[it <a\<y s ] 



Is 1 ls2 
Is: ) ) • • • 

7s 7*i 



(28) 



since, according to the result to be established below, isi/ls, ls 2 /ls 1 , ■ ■ ■ are independent of 
7 S and 74. It follows that the gamma process is Markovian. A similar argument shows that 
the gamma bridge has the Markov property. In particular, we have 



lt_ 

It 



< a 



Is Isi ls2 

It' It' It 



lt_ 


< a 




_lT 




~lt_ 


< a 




.It 





Is ls\ ls2 

It Is' 1 S1 ' 
7s 
It. 



(29) 



since the random variables i si /ls, ls 2 /lsi, ■ ■ ■ are independent of 7i/7T and i s /it- 
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Proposition 2 Let { 7 t}o<t<oo be a standard gamma process. Then for T > t > the 
random variables jt/lr an d 7t ^re independent. 



Proof. For the joint distribution of these random variables let us write 



It 



F(y,z) = 

We note that this can be rearranged in the form 
F(y,z) = 



— < y n 7t < z 

It 



(30) 



y 



It < (7t -It) n 7* < z - (7t - 7*) 

1 - y 



(31) 



Conditioning with respect to 7r — 7t, we use the fact that 7 t and 7 t — 7t are independent, 
and that j t has a gamma distribution with parameter mt, to deduce that 



F(y,z) 



T[mt] 



E 



(32) 



Differentiating each side of this relation with respect to y and z, we obtain the following 
expression for the joint density function: 



f(y,z) = 



l 



T[mt] 



E 



Jo (1 



It ~lt 



yf 



5 x 



y 



i-y 



(7t - It) 



x % - [z - ( 7T - 7i )])a; m *- 1 e- :E dx 



(33) 



Here we have used the relation d x \s x < a \ = —5(x — a), where {£(z)}zeR denotes the Dirac 
distribution. Integrating out the first delta function we thus have 



f(y,z) 



T[mt] 



E 



It - It 



[i-y) 2 \i-y 



ilT-lt) ~\z- (7r-7t)] 



(7t - It) 



mt—l 



,"i(7T-7t) 



(34) 



for < y < 1 and z > 0; and hence after some rearrangement we obtain 



f(v,z) = 



jmt— l^j y^—mt—1 



T[mt] 



E 



(7t - 7t ) m V^ (7T ~ 7t) 5( 



i-y 



(35) 



Now we introduce the Fourier representation 

5(x) = — / e iXx d\ 

for the delta function, interpreted in a distributional sense, from which we deduce that 



(36) 



f(y,z) 



yint-ln _ y\-mt-l -y 



F[mt] 



2tt 



(7t - 7t ) m 'e"^ (7T - 7t) e iA ^ (7T - 7 * ) ] dA. (37) 



10 



Writing E for the result of the expectation appearing inside the integral above, and making 
use of the fact that 77- — j t is gamma distributed with parameter m(T — t), we have 

1 f°° 
T[m(T-t)]J 



1 



00 



T[m(T-t)] J 

(i-y) mT 



T[m(T-t)\ J 
{l-y) mT T[mT] 



u mT-l e -u e iXu du 

(38) 



T[m(T-t)] (l-iA) mT ' 

where we have made use of (|7J) to deduce that the characteristic function of 7^ is 1/(1— iA) mT . 
Substituting (1381) into (1371) we obtain 

'<**> = r M r!l(T- t) ] - (T^p e-dA, (39) 

and hence 

„.mt— If-] „.\m{T—t)-l r tnT-l p -z 

f(y, z) = y % ^ V^r ( 4 °) 

i5[mt,m(l —t)\ 1 [mi J 

for < y < 1 and z > 0. Thus we see that the joint density for 7j/7t and jt factorises into 
the product of a beta density for 7j/7t and a gamma density for j T , as desired. □ 

The result of Proposition [2] is a special case of the following more general result: 

Proposition 3 Let {7t}o<t<oo be a standard gamma process. Then for T > u > t > the 
random variables (j u — 7*)/(7t — It) an d It — It o,re independent. 

Clearly, Proposition [2] follows as a special case of Proposition [31 The following lemma is 
a classical result (Lukacs 1955, Yeo & Milne 1991) which can be used as the basis of a proof 
of Proposition [31 

Lemma 1 Let A and B be independent gamma- distributed random variables with parame- 
ters p and q, respectively. Then A/ (A + B) and A + B are independent, A/(A + B) has a 
beta(p, q) distribution, and A + B has a gamma(p + q) distribution. 

Proof. For independence it suffices to show that the joint Laplace transform of A/ (A + B) 
and A + B factorises. In particular, for positive a, /3 we have 

poo poo p-\ -a ig-l -b 

L J Jo Jo r[p] T[q] V > 

Setting x = a/(a + b) and y = a + b, we have a = xy and b = (1 — x)y, and hence 
da db = y dx dy. We see that 

pi ~p-iyi _ ™\<?-l p°° „.p+q-l„-y 

1 J Jo B[p,q] J T[p + q] 

= E [ e -^/(A+B)^ E [ e -/?(A+B)] (42) 
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It follows that A/(A + B) and A + B are independent and have the distributions stated. □ 

The proof of Proposition [3] follows if we set A = 7„ — j t and B = jt — Ju- A proof of 
Proposition [2] is obtained if we set A = j t and B = jt — It- 



V. VALUATION OF AGGREGATE CLAIMS 

Our objective is to calculate the value at t of a contract that pays Xt at T. We assume 
that Xt is strictly positive and integrable. For simplicity of exposition, in this section 
we take Xt to be a continuous random variable; the adjustments required for the more 
general situation are straightforward. We assume that the default-free interest rate system is 
deterministic, that Q is the risk-neutral measure, and that the market filtration is generated 
by an aggregate claims process {£t}o<t<T of the form £ t = XTjtT, where {jtr} is a standard 
gamma bridge under Q, with parameter m, which we take to be independent of Xt- The 
value St of the contract at t < T is given by St = PtyE [Xt \ J 7 t), where Tt — o ({£s}o<s<*)- 

Proposition 4 The aggregate claims process {£t}o<t<T has the Markov property. 

Proof. For the Markov property we must verify that 

Q [& < a | F s ] = Q [£ t < a | £ s ] (43) 

for all s, t such that < s < t < T. It suffices to establish that 

Q & < a | ^, U, • • • > U) = Q[£t<a\ 6] (44) 

for all t > s > Si > S2 > ■ ■ ■ > s n , and for all n > 1. We use the representation 
{7tr} — {it/lr}, where {74} is a standard gamma process with rate m. Then we have 



[it < a I 



S21 



£t<a 
it<a 



v 7« y 7«1 7s 2 
, J\-T j J\-T • 



X, 



It It 
Is 7*1 7s 2 



7r 



(45) 



7t 7s 7sx 

But 7ai/7s) 7s 2 /7sd • • • are independent of ^ and £ s , which gives us the desired result. □ 

By virtue of the fact that {^} has the Markov property and that Xt is jF T -measurable 
we are able to simplify the expression for St so that it takes the form 



S t = PtTE[X T \& 



(46) 



The conditional expectation appearing here can be carried out in closed form, leading to the 
following pricing formula: 

Proposition 5 The value St at time t < T of the aggregate claim that pays the continuous 
random variable Xt > at time T is given by 



S, = P, 



x)x 2 ~ mT (x - ^ t ) m ( T - t )- 1 dx 
tT j^p{x)x 1 - mT {x - ^)™( T -*)- 1 dx' 



(47) 



where {p(x)} 0<x<oc is the probability density of X T . 
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Proof. The conditional expectation (1461) can be written in the form 

POD 

E [X T | £t] = / x-n t (x)dx, 
Jo 



(48) 



where {ir t (x)} is the conditional density process for Xt, which by virtue of the Markov 
property of is given by 



7r t (x) = ^-Q[X T <x\Ct}. 



(49) 



We can compute Tr t (x) by use of the following form of the Bayes formula: 

^ X) ~ j-p{x)p&\X T = xW { ) 

where p (^ | Xt = x) is the conditional density for £ t , valued at Pf Specifically, we have 

p(t\X T = x) = ^-®[Zt<Z\X T = x] 



d£ 
_d 



[x T j tT < £ | x T = x] 



7tr < - 

X 



(51) 



Therefore, writing {f(y)}o<y<i for the density function of the random variable jtr we find 



d r /x 

p(i\X T = x) = - / f(y)dy 
d£ Jo 



-/(-)■ 

X \x 



(52) 



Hence by Proposition [T] we have 



n(P\X r\ S (^) mt - 1 (l-e/a:) m(T - f) - 1 



i{x>£}£ 



mt—l 



X 



l-mT/ x _ rtm(T-t)-l 



B [mi, m(T - *)] 

The conditional probability density function {7Tt(a;)}o<t<T,x>o for -^x is thus given by 

P (x)x 1 - mT (x - ^y^)- 1 
TT t {x) - 1 {X> ^ } j°o p ^ xl - mT ( x _ 6)m( T-t)-l dx ' 

from which the desired result (|47l) follows at once. 



(53) 



(54) 
□ 
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With these results at hand we are also in a position to price a simple stop-loss reinsurance 
policy. For such a policy the value process is given by (2), and hence we have 



poo 

C tT = PtT / (x - K) + ir t (x)dx 
Jo 



f™(x - K) + p{x)x 1 - mT (x - it) m{T - t] - 1 dx 
= PtT J^p(x)x 1 - mT (x - dx ' (55) 

It should be evident that once a time t has been reached such that £ t > K, then C u t = 
PuT(St — K) for all u such that t < u < T. In other words, once a sufficient number of 
claims have accumulated the option is sure to expire in-the-money. 



VI. VALUATION OF GENERAL REINSURANCE CONTRACTS 

In the previous section we showed how one works out the reserve process for an aggregate 
claim that pays Xt at T, and we were also able to determine the value process of a stop-loss 
contract that pays (Xt — K) + at T. In this section we consider the more general situation 
of a contract that at a fixed time t < T allows the policy holder the option of commuting 
the claim Xt in exchange for a pre-fixed settlement K. Let us write Co t for the value at 
time of such an option; then clearly we have 

C ot = P ot E [(S t - K)+] , (56) 

where St is the value at t of the claim that pays Xt at T. With reference to Proposition [5l 
it will be useful to introduce a function S(t, y) for < t < T and y > by setting 

f™p(x)x 2 - mT (x - yY^-^dx 
V) = PtT p(x)xi~™ T (x - y)m(T-t)-i dx ( 5? ) 

Then the value of the claim is given by S t = S(t,£ t ), and the value of the option can be 
written in the form 

C t = P tE[(S(t,Zt)-K) + }. (58) 

Since the payout of the option is a function of £ t , one way of working out the expectation in 
( Jo's"]) is to obtain an expression for the price A ot (y) of an Arrow-Debreu security that pays 
5(£ t — y) & t t, where y > is a parameter. Thus we have 

A ot (y)=P ot E[d(Z t -y)}, (59) 

and for the option we can write 

POO 

C ot = / A ot (y)[S(t,y)-K] + dy. (60) 
Jo 



We shall calculate A ot (y) and use the result to determine the expectation (I5"8"|) . We state 
the result of this calculation first, the proof of which is given at the end of this section. 
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Proposition 6 The price A ot (y) at time of an Arrow-Debrue security that pays <5(£ t — y) 
at t is given by 



A ot (y) = ^ 



y 



mt—l 



Of 



B[mt,m(T -t)) 



p(x) x^^ix-yr^-^dx. 



(61) 



By comparing (l57j) and (I6T1) we observe that the integral term in (1611) cancels with the 
denominator in the expression for S(t,y). After some rearrangement we thus obtain 



G 



Of 



Poty 



mt — 1 



B[mt,m(T -t)} 



p{x) (xP tT - K) x x - mT {x - y) m(T -* ) - 1 dx 



<- J y 



dy. (62) 



We are now left with the task of finding the critical values at which the argument of the 
max-function in the integrand of ( 1621 vanishes. Suppose that S(t, y) is monotonic in y; then 
there is at most a single critical value y*, obtained by solving the following equation: 



p(x) [xP tT - K) x x - mT {x - t/*) m(T_t)_1 ^ = 0. 



(63) 



The lower limit of the outer integration in the expression for Cot above can then be changed, 
and we have 



Of 



Poty 



mt—l 



p(x) {xP tT - K) x^ mT (x - yr^-'dx 



r B[mt,m(T -t)] 

This expression simplifies further if we swap the order of integration as follows: 
Pot 



dy. (64) 



Of 



B[mb,m(T-t)] 
Making the substitution y = xz we then obtain 

Pot 



poo px 

/ / p{x) (xP tT - K) y mt - x x l ~ mT [x - y^-^dy dx. (65) 

Jy* Jy* 



Of 



B[mt,m(T -t)] 



/ p{x) {xPrr - K) / z mt ~ l {l - z^'^dzdx. (66) 

Jy* Jy*/x 



Let us now introduce the complementary beta distribution function B(u) with parameters 
mt and m(T — t) by the expression: 



B{u) 



f l uZ mt-l( l _ z yn{T-t)-l dz 

J 1 ^*- x (l- z) m ( T -Q- l dz 



(67) 



We call this the "complementary" distribution because the integration ranges from u to 1. 
Clearly, the denominator in fl67|) is B[mt, m(T — t)]. We thus find that the integration over 
the variable z in fl66|) combines with the factor B[mt, m(T — t)] appearing of that expression 
to give a cumulative beta distribution function, and for the option price we have 



POO 

Cot = Pot / p(x) (x P tT - K) B(y*/x) dx. 

Jv* 



(68) 
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We remark, incidentally, that a sufficient condition for S(t, y) to be monotonic in y for 
fixed t is m{T — t) > 1. To see this, we differentiate S(t,y) with respect to y, assuming the 
stated condition, and after some rearrangement we obtain 



dS(t,y) 
dy 



( 



P tT [m{T - t) - 1] 



J°° p(x)a 2 (x)dx J°° p(x)(3 2 (x)dx 



\ 



V 



(/J* ' p(x)a(x)P(x)dx 



- 1 



(69) 



/ 



where a 2 (x) = x 1 mT ( x — y) m ( T *) and /3 2 (x) 



x 



1—mT 



{x- y ynp-t)-* m \fm{T-t) > 1, then 



the integrals exist, and it follows on account of the Schwartz inequality that dS(t, y)/dy > 0. 



Proof of Proposition O It suffices to determine the expectation E[5(£ 4 — y)\. By use of 
the Fourier representation (1361) we can write 



(70) 



E[*(£t - V)] = =- / e^E[e iA ^]dA. 

^ J—oo 

Since Xt and 'jtr are independent, it follows from the tower property that 

E[e iA6 ] = E[E [e iAXT7tT |X T ] 

p(x) E [e iAx ' 7 ' T ] dx 

3 

p(x) 4>(\x) dx, 

where <p(v) = E [e li/7tT ] is the characteristic function of 7^. We deduce that 



(71) 



1 poo poo 

-y)} = 7r e ~ iXy / 

/7r J -00 Jx=0 



p(x) 4>(Xx) dx d\. 



(72) 



Thus, by interchanging the order of integration and using the fact that the inverse Fourier 
transform of the characteristic function is the density function we have 



p(x) 



x=0 

00 



1 

2^ 



e- iXy <P(\x) d\ 



dx 



p(x) — 

x=0 x 

00 



/ Q-ivy/x, 

2tt 



M dv 



dx 



I p(x)~f (-) dx, 
Jo x \xJ 



(73) 



where / is the density function of j t T- Substituting the expression (jTTj) for / into (J73]) we 
find that 



r°° 1 
E[5(Ct-y)] = / p{x)-t {x>y} 

Jo x 



(y/x)™*- 1 ^ - y/x)" 1 ^-^- 1 



dx 



,mt—l 



B[mt,m(T -£)] 



B[mt,m(T -t)] 
p{x) x l - mT (x-y) m P-V- l dx, 



(74) 
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which verifies the claim. □ 
We remark that the price of the Arrow-Debrue security can be put in the form 

Cp{y/u)u mt - 2 (i - u) m ( r -*)- 1 d'u , N 

Aot(y) = Pot ° 1 , 75 



by use of which the normalisation J °° Aot(?/)d?/ = Pot can be checked. It follows also from f!75|) 
that the characteristic function $^(A) of £ t is given by the beta average of the characteristic 
function <fi x of X T : 

VII. DISCRETE CASH FLOWS 

Thus far we have considered the case for which the terminal cash flow is a continuous 
random variable. In this section we consider the example for which Xt takes values in 
a discrete set {:£i}i=i,...,n- The corresponding a priori probabilities will be denoted {pi}- 
The calculation presented in Section |V] holds and we obtain, instead of (HTj) . the following 
expression for the value process: 

q _ P Y,iPiX 2 i~ mT (xi - 6) m(T " f) ~ 1 i{6<^} (77 , 

It is straightforward to verify that expression (ITT]) converges to the correct terminal value 
as t approaches T. To see this, suppose that for some u G Q the value of Xt is Then 
for that choice of u we have 

and hence, after some rearrangement, 

/ \ l-m(T-t) 

~ ~t 1-mT ( l- 7tT y- m ( T -*) ' ™ 

It follows at once that St = Xk- 

We proceed now to value a reinsurance contract that pays (St — K) + at time t. For this 
purpose we need the price of an Arrow-Debreu security with payoff 5(£t — y) at t. In the 
discrete case the Arrow-Debreu price is given by 



yTTlt—l 



(80) 



Substituting flHUl) and the function 



s & y) = p ^^\- mT r " (r t) ,7 l>y} ( 81 ) 
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FIG. 1: The value function S(t,y) associated with the reserve price in the case of a discrete cash 
flow at time T taking four possible values. The parameters are chosen such that {x±, X2, X3, X4} = 
{1,2,3,4}, {pi,P2,P3,PA} = {0.5,0.2,0.2,0.1}, m = 2.0, r = 5%, and T = 1. For a given time t 
the value function represents the reserve required if the aggregate claims amount to y. 



into (160]) we obtain, after some rearrangement, 



Of 



Of 



B[mt,m(T -t)] 



.mt—l 



8=1 



Pi xl~ mT ( Xi - yr^-H^yyiPtrX, - K) 



dy. (82) 



A discrete version of formula ( |69l) shows that 5(t, is increasing in y if m(T — t) > 1, and 
decreasing in y for y G (xfe, x^+i) for each fc = 1, . . . , n — 1 if m(T — t) < 1. See Figured] for 
the typical behaviour of S(t, y) when Xt takes four possible values. For fixed t there is at 
most a single critical value y = y* for which S(t, y) = K, when y 7^ xj~ for all fc. We thus 
have three scenarios to consider, namely: (I) S(t,y) is increasing in y at y — y*\ (II) the 
critical value y* is at y = Xk for some k; and (III) S(t,y) is decreasing in y at y = y*. 

We therefore analyse the price of the reinsurance contract in these different scenarios. In 
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case (I) the integrand in (I82I) is nonzero when y G (y*, oo), and we have 

n 00 

Cot = B[mt ^[ T _ t) { Ev^- mT {PtTXi ~K)j y^\ Xi - y)^l { ^ ( 83 ) 

«— ^* 

The ?/ integration in fl83|) can be carried out by observing that for X{ > y* we have 

r Xi f 1, \ mt ~ l f i,\ m ( T -*)- 1 
y^^-yr^-H^ydy = xT T ~ 2 / ( -M 1 - I) dy 



X{ J \ Xi 

mT—1 I „mt—l 



f z mt -\i - z^-^dz, (84) 

Jy*/xi 



where we have made the substitution y = X{Z. Therefore, the price of the reinsurance 
contract can be expressed in terms of the complementary beta distribution function with 
parameters mt and m(T — t): 

n 

l>y *}Pi(PtTXi - K)B(y*/xi). (85) 

i=l 

If there is no critical value in the range (xk,Xk+i), then y* = Xk for some k. Hence the 
pricing formula in case (II) is identical to the result obtained in ( |85l) . with y* = Xk- In 
case (III) there are two distinct regions for which the integrand in fl82|) is nonzero. These 
are given by y G [xk,y*) and y > Xk+i for some k. Hence the pricing formula is similar to 
that obtained in fl85l) . except there are additional terms arising from the integration over 
the range [x^, y*). 

As an example of a discrete cash flow we consider the binary case where Xt can take the 
values x , x\. In this situation the critical value y* < x can be worked out by solving 

Po(PtTX - K)xl mT (x - yT 1 ^' 1 = Pi(K - P tTXl )x\~ mT {x x - y*) m ( T -^ (86) 
for y*. A short calculation shows that 



0x\ - x 

— , where 6 

9—1 



Px(K - PtTXl) ( Xx 



1-mT 



m(T-t)-l 

(87) 



Po(PtTX - K) \x 

It follows that the price of a reinsurance contract in the case of a binary payoff is given by 
Cot = Po(Potx - PotK)B(y*/x ) + Pi(P rXi - PotK)B{y* / x x ) . (88) 

VIII. OPTION PRICE PROCESS 

We generalise now the analysis of Section [VTl to derive an expression for the price process 
of a call option on the value of the reserve St at time t associated with the claim Xt- As 
before, we let K be the strike. Then the value of the option at time s < t is given by 

C st = P st E[(S(t,£ t )-K) + m. (89) 
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Once again we find it convenient to obtain first the price process for the Arrow-Debreu 
security. This is on account of the relation 



C 



St 



IJo 

oo 



6 



Sfa-y) (S(t,y)-K) + dy 

POO 

P st E[5(^-y)m(S(t,y)-K) + dy 
Jo 

A st (y) (S(t,y)-K) + dy, 



where S(t,y) is defined as in ( 1571) . and {A st } < s < t <T is given by 

A st (y)=P st E[5^ t -y)\Q. 
By taking the conditional expectation we obtain the following result: 



(90) 



(91) 



Proposition 7 The price process {A st (y)} < s < t < T of the Arrow-Debreu security that pays 
out 5(£ t — y) a t t is given by 



A st {y) = P t 



; B[m(t - s), m(T - t)} /~ p(x) x l ~ mT (x - Q^-^dx' 
where y > £ s and {p(x)} is the probability density of Xt- 



(92) 



This result is established later in this section. By substitution of fl92l) in (|90|) we see that 
the price process of the option is given by 



a 



St 



St 



B[m(t - s), m(T - t)} /°° p(x) x 1 ~ mT {x - Q^-^dx 



X / (y - ^(t-s)-l 



p{x)x 1 - mT {x - y^P-^iPtTX - K)dx 



dy. (93) 



Assuming that there is only one critical value y* that solves fl63l) . we find that the integration 
over y in (193]) vanishes for y smaller than y*. In this case, we can lift the max-function in 
the integrand, and by interchanging the order of integration we obtain 



C. 



P. 



si 



st 



B[m(t - s),m{T - t)) L°° p(x) x l ~ mT (x - Q m( ?- s )- l dx 



xP 



tr 



x=y~ 



p{x)x l ~ mT {P tT x -K) (y- ^-^{x - y)^-*)- 1 dy dr. (94) 



y=y 



Let us analyse the y integration. Making the substitution y = z + £ s we find that 



_ p ym{t-s)-\l x _ y \m(T-t)-l 



(y - C) 



dy 



'y=y 



z 



m{t—s)—l 



X 



- C s - z) m ^- l dz. (95) 



«=2/-6 
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A further change of variable obtained by setting w = z/(x — £ s ) gives 
' 8 z m ^-\x - & - z^-^dz 



x _Qm(T-.)-l f w^-^il-W^-^dw. (96) 

J w= y^is 



We see that together with the beta function in the denominator of (1941) the integral term in 
the right side of ( 1961) gives rise to a complementary beta distribution function. Therefore, 
the call price can be written in the form 



Finally, we observe that the quotient in the integrand is the conditional density tc s (x). The 
call price at time s < t thus reduces to the following expression: 

C st = P st n s (x)(xP tT — K)B ( J dx. (98) 

J x=y* V Ss / 

As in the case of the initial price of the option, the range of integration in (I9"8"j) must be 
modified appropriately if there is more than one critical value for which f)63p is satisfied. We 
now proceed to derive the expression for the Arrow-Debreu price process. 

Proof of Proposition By use of the Fourier representation (1361) we have 

1 f°° 

E im - y) \ts] = ^J e- iX y E [e iA ^ | &] dA. (99) 

To determine the conditional expectation Ee lA ^|^ s ] the following result is handy: 

Lemma 2 Let {£,t}o<t<T be given by £ t = XxjtT, where {jtr} is a gamma bridge and Xt is 
an independent positive random variable. Then for fixed s such that 0<s<t<Twe have 

6 = 6 + Z T 5 tT , (100) 
where Zt = (1 — J s t)Xt, and where the process {5tr} s<t<T , defined by 



ht 



101) 



1 - IsT 

is a gamma bridge over the interval t e [s,T] and is independent of 6 o,nd Z?. 

By use of fllOOp and the tower property we find that 

E [e iA?t | 6] = E [ e iA(fs+ZT5tT) | 6] 
= e iA6 E [e iAZr<5tT | £j 

= e iA ^E [E [e iXZ ^ |&, 5 tT ] | £ s ] . (102) 



21 



Since Z T = X T — £ s , and since {5 t r} is independent of £ s and X T , the inner expectation can 
be carried out explicitly by use of the conditional density for X?, and we obtain 



" POO 

E[e iA€t |^] = e iA « s E / e^-^^vr.^dx 

lJx=£ s 



,Jx=£ a 

By substituting (11031) in fl99l) we deduce that 



(103) 



1 /*oo /*oo 

E [5(6 - y)|&] = 7T / e" iA ^) / $ 5 [A(x - &)] vr,(x) dxdA 

/•CO / -i /*GO \ 

= J e- iX{y -^MHx-Zs)]d\jir s (x)dx, (104) 

where $5 is the characteristic function for Str- By use of the substitution z = X(x — £ s ) we 
then find that 

E[*(6-y)|&] = f° fe-^^^^-Tr^di 

*- (a0 -AWW (105) 



00 



£s \% S. 



where /<5 is the probability density of 5 t T- Since <5 4 t is beta distributed with parameters 
m(t — s) and m(T — t), we deduce, after some rearrangement, the expression obtained in 
( |92|) for the Arrow-Debreu price. □ 

Proof of Lemma [2j The decomposition (llOOp can be verified by direct calculation if one 

sets {-ftr} = {it/ It}, where {74} is a standard gamma process. To see that {5 t T}s<t<T is, 
for fixed s, a gamma bridge over [s,T] it suffices to note that S t r = (it — J a )/(jT ~ 7s) and 
that {74 — 7s} s <i<oo is a gamma process. In particular, we observe that the independent 
increments property holds, and that jt — 7s is gamma distributed with mean m(t — s). 
Finally, to see that {5tr} is independent of £ s and Zt it suffices to show that 5 t T, 7s and 
7t — 7s are independent. We have: 

Q({5 tT < a} n {7t - 7s < b} n {7s < c}) = E[l {(5tT<a} l {7T _ 7s<b} l {7s<c} ] 

= E[E[l {5tT<a} l {7T _ 7s<6} l {7s<c} I 7 S ]] 

= E[E[l {5tT<a} l {7T „ 7s<f)} I7J l{ 7s<c }] 

= E[E[l {(5tT<a} l {7T „ 7s<f)} ] l{ 7s<c }] 

= E[l {(5tT<a} l {7T _ 7s<6 }] E[l {7s<c} ] 

= E[l{,5 tT<a }] E[l{ 7T _ 7s<b } ] E[l {7s<c} ] (106) 

In going from the fourth to the fifth line we have used the fact that 7 S is independent of 5 t T 
and 7t — 7s; which can be checked by use of the independent increments property of {7*}. 
In going from the sixth to seventh line we have used Lemma [1] together with the fact that 
we can write 8 t r = B/(A + B) and 7t — 7 S = A + B, with A = 7^ — j t and B = j t — j s , 
from which it follows that 5 t r and 7^ — 7s are independent. □ 



The result of Lemma [2] leads to the following observation concerning the model calibra- 
tion. Suppose that the aggregate claims process is given, and that we reinitialise the model at 
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some specified intermediate time. We would like the dynamics of the model moving forward 
from that intermediate time to be consistently represented by an aggregate claims process 
of the same type. Indeed, it follows from Lemma [2] that the process {r] t } s <t<T defined by 

7} t = Z T 5 tT (107) 

is an aggregate claims process spanning the time interval [s,T]. The random variable Zt 
can be thought of as representing the information about X? that is "not yet revealed" at 
time s. The idea is that at time s the value of £ s is known, and the "new" gains process 
{rjt} s <t<T begins to reveal the value of Z T in such a way that r/ s = and r\ T = Z T . 

Alternatively, at time s we can use the knowledge of £ s to compute the "new" a priori 
density for Xt- Thus, at time s the a priori density p(x) for Xt is replaced by the appropriate 
a posteriori density tc s (x). On account of the relation Zt = Xt — C, s we have 

Q[Z T < *|&] = Q[X T < z + £ s |£ s ], (108) 

from which it follows that the conditional density of Zt is given at time s by tt s (z + £ s ). 
We can think of 7r s (£ s + z) as a "new" a priori density, now for the random variable Zt- 
Given this density we calculate the conditional probability Q[Z T < z\rj t ] for t G [s,T]. By 
the method used to establish Proposition H] and the probability law for the gamma bridge 
{5tT} we deduce that the associated density function is given by 

from which we see that the value process can be represented in the following form: 



S t 



tT 



(110) 



Making the substitution z = x — £ s and also substituting r] t = £ t — this expression reduces 
to the value process obtained in (H7|) . 



IX. EXAMPLE: GAMMA-DISTRIBUTED CASH FLOW 

When the terminal payout X T of the cumulative gains process (j3J) is gamma distributed 
with mean nmT and variance K 2 mT for some choice of k, the resulting value process {St} 
has an especially simple structure. In particular, we are lead back to the "Q-gamma" model 
discussed in the introduction. This can be seen as follows. Let {7t} be a standard gamma 
process with rate m, and let {7^} be the associated gamma bridge. Then Xt and k^t have 
the same distribution; but since jt and {jtr} are independent, it follows that {X T 7 iT } and 
{/t7T7tr} have the same probability law; therefore, {£t} has the same law as {^74}, and 
hence is a Q-gamma process, with scale k and standard growth rate m. The fact that £t 
is gamma distributed can be verified directly as follows. The characteristic function of Xt 
is cj)x(X) = (1 — in\)~ mT . Substituting this into (1761) and setting z — (1 — u)/(l — in\u), 
we deduce that $^(A) = (1 — i«;A) _m *, which is the characteristic function of a gamma 
distributed random variable with mean nmt and variance K 2 mt. 
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It is interesting to note that although the Q-gamma process has independent increments, 
the cumulative gains process ([3]) has dependent increments. In particular, for the covariance 
of £ s and £ t — £ s in the general case we have 



Cov[&,&-&] 



T(mT + 1) 



s(t - s) 

J>2 



(E[X T ]Y 



111) 



Hence a necessary condition for independent increments is given by (E[Xt]) 2 = mT VarLYr]. 

We conclude the paper by working out in some detail the value processes for various 
claims in the Q-gamma model. For the density of X T we have gr T (x), where gr t ( x ) is 
defined in ( 1T4|) . Substituting the expression for the density function into ( H71) and carrying 
out the relevant integration, we are led to the following expression for the reserve process: 



S t = P tT (Z t + Km(T-t)). 



;ii2) 



Therefore, {S t } in this case is a linear function of {£*}■ We observe that So = Pqt^^T and 
that St = Xt, as required. An alternative derivation of (11121) is as follows. Since {^} is 
a gamma process with scale parameter k and standardised growth rate m, by the Markov 
property we have St = PtT^>[(,T\(,t], and ( 11121) follows immediately as a consequence of the 
independent increments property of the gamma process. 

These relations lead to simplifications in the valuation of contingent claims. Let us work 
out, for example, the value Ctr at time t of a simple stop-loss reinsurance contract that pays 
out max(Xr — K, 0) at T for some fixed threshold K. In the Q-gamma model we have 



C tT = P tT E[^ T -K) + \a 
and hence by use of the independent increments property we deduce that 

poo m(T-t)-l -z 

C tT = P tT / (kz + it- K)=^-= — dz 

T V\m.(T 

= PtT 



113) 



(K-k)/ K T[m(T-t)] 



J[m(T-t) + l,(K-£ t )/ K ] r[m(T-t),(K-Z t )/ K ] 
« — — — ~ (A - — 



T[m(T - t)] 



T[m(T-t)] 



114) 



where F[a, z] = x a l e x dx denotes the incomplete gamma integral. 

We proceed to calculate the associated Arrow-Debreu price A st in this model. By substi- 
tuting ffTT2"|) in ([22]) we deduce that 



-rn(t-s) / i 



(115) 



It follows by use of flll2p that the price at time s of a reinsurance contract with payout 
{S t - K) + at t is 

C st = P st E s [(S t - K) + ] 

POD 

= / A st (y)[P tT (y + Km(T~t))-K] + dy 
Jo 

T[m(t - s) + 1, K- 1 ^} _j T[m(t - s), n~ l R s 



p. 



sT 



where R q 



T[m(t - s)) 
P~ T X K - (S s + Km(t- s)). 



T[m(t - s)\ 



(116) 
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